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Abstract. It is known that the class Up, of generalized s-selfdecom- 
posable probability distributions, can be viewed as an image via 
random integral mapping of the class ID of all infinitely di- 
visible measures. We prove that a composition of the mappings 
jPx ^ j p 2 ^ _ ^ jPn j s a g a i n random integral mapping but with a 
new inner time. In a proof some form of Lagrange interpolation 
formula is needed. Moreover, some elementary formulas concern- 
ing the distributions of products of powers of independent uni- 
formly distributed random variables as established as well. 
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Abbreviated title: On generalized s-selfdecomposability 

Let a real /3 and infinitely divisible probability measures Uj (on Euclidean 
space lR d or a Banach space E) be such that 

Ti/ n (vi * u 2 * ... * v n )* n P * o~x n => A* , as n->oo. (*) 

Then we write \i G Up and call \i a generalized s- self decomposable measure. 
In a series of papers Jurek (1988, 1989), Jurek and Schreiber (1992) it was 
proved, among others, that for not degenerate \i in (*) we must have (3 > — 2 
and that Up (/3 > —2) form an increasing family of convolution semigroups 
that " almost" exhaust the whole class ID of all infinitely divisible measures. 
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However, for the purpose of this paper the most crucial is the fact that gener- 
alized s-selfdecomposable measures admit the random integral representation 
(1), (see below), which is a particular case of the following representations 



where I = (a, b] C IR + , h : IR + — > R, Y u (-) is a Levy process and, 
r : M + — > M. + is a monotone function (deterministic time change in Y u ) . 

In fact, it was proved that many classes of limit laws can be described as 
collections of probability distributions of random integrals of the form (**). 
Later on, this lead to the conjecture that all classes of limit laws should 
admit random integral representation; cf. Jurek (1985; 1988) and see the 
Conjecture on www.math.uni.wroc.pl/~zjjurek0. More recently the method 
of random integral representation was used among others by Aoyama and 
Maejima (2006), Maejima and Sato (2009). Here we give yet another example 
of calculus on some infinitely divisible laws. 

In this paper we will prove that the class of the integral mappings (1), 
for (3 > 0, is closed under compositions, that is, their compositions are of the 
form (**) with the properly chosen time change r; (Theorem 1, Proposition 
2). As an auxiliary result we found a decomposition of number 1 as a sum 
of products of complex fractions; (Proposition 1). Also compositions of the 
mappings (1) are described in terms of the Levy-Khintchine triples: (The- 
orem 2). Auxiliary Lemmas 2 and 3 give probability distribution functions 
(p.d.f.) of products of powers of independent uniformly distributed random 
variables as linear combinations of other p.d.f. 

1. Introduction and main results. 

The results here are given for random vectors in M d . However, proofs 
are such that they are valid for infinite dimensional separable Banach spaces 
E when ones replaces a scalar product by the bilinear form on the product 
space E x E, where E denotes the dual space. Of course, (M. d ) = R d . 

Throughout the paper C(X) will denote the probability distribution of 
an M. d - valued random vector X; (or a real separable Banach space E- valued 

1 It might be of an interest to recall here that S. D. Chatterji's subsequence princi- 
ple claiming that: Given a limit theorem for independent identically distributed random 
variables under certain moment conditions, there exists an analogous theorem such that 
an arbitrary- dependent sequence (under the same moment conditions) always contains a 
subsequence satisfying this analogous theorem was proved by David J. Aldous (1977). Al- 
though, we do not expect that the above Conjecture and Chatterji's subsequence principle 
are mathematically related, however, one may see a "philosophical" relation between those 
two. 
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random element X if the Reader is interested in that generality). Similarly, 
by Y u (t), t > 0, we will denote an IR d - valued (or an E- valued) Levy stochastic 
process such that C(Y V (1)) = v. Recall that by a Levy stochastic process 
we mean a process with stationary independent increments, starting from 
zero, and with paths that are continuous from the right and with left limits 
(that is, cadlag paths). Of course, v e ID, where ID stands for all infinitely 
divisible measures on M. d (or on a Banach space E). 

For j3 > and a Levy process Y v (t),t > 0, we define mappings 

J^\v) = j\v): = C{\ tW dY v (t))=£{[ tdY v (t P )) (1) 

■/(o,i] Ao,i] 

and the classes Up : = J 13 (ID). To the distributions from Up we refer to as 
generalized s- self decomposable distributions. 

These classes of probability measures were originally defined as limiting 
distributions in some schemes of summations; cf. Jurek (1988 and 1989). In 
particular, the class U = U± of s-selfdecomposable was defined by non-linear 
shrinking operations U r ,r > 0, (for x > 0, U r (x) := max(0,x — r)); cf. Jurek 
(1981). 

Remark 1. Since the process Y has values in a metric separable complete 
space we may and do assume that the paths of Y are cadlag; cf. Theorem 
A. 1.1 in Jurek and Vervaat (1983), p. 260. Since the random integral in (1) 
is defined by a formal integration by parts formula, therefore the random 
integral in question do exist; cf. Jurek- Vervaat (1983), Lemma 1.1 or Jurek- 
Mason (1993), Lemma 3.6.4, on p. 119. 

For a positive natural m and a sequence of positive real j3i, 02, Pm and 
a probability measure v G ID, let us define the mappings 

;= jP m ( J Uh,...An-i}( J ,)) =C ([ t dY j{01 ,..., 0n _ l} (t^)) . 

J (OA] 

Our main results say that the above compositions can be written as an 
integral of the form (**) with a suitable chosen time change r. Furthermore, 
compositions are expressed in terms of the individual random integrals. 

Theorem 1. For positive integers f3\, f3 2 , fim an d an infinitely divisible 
probability measure v we have 

jVu~M {v) =C {! t dY v (r { p u ...,p m} (t)) = ll^- M {v) (2) 
.7(0,1] 
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and the time scale change r{p u _$ m } is given by 

r^,..., M (t) = F[Ul /h ■ U 1 ^ • ... • UU^ < t], < t < 1 , (3) 

where Ui 's are stochastically independent random variables uniformly dis- 
tributed over the unit interval. 

If all j3i, (3 n are different then 

n n n 

r{/W»>(*):=£ C *»* ft ' C itn {fi u ...^ n ) = C itn := \\ -=-^-=- (4) 

and, in particular, we get the equality: YTj=i Cj,n = 1 ■ 
If Pi — P2 — ■■■ — fim = ot ( m > 1) then 

r { a,...,a } (t)=t a Y, ° g . /or < t < 1 , (5) 



m — times 



3=0 



(Other collections of /3's are discussed below.) 

Remark 2. Note that if £(A) denotes the exponential random variable with 

the parameter A then e - ^^ = U l ^ x . Hence if £i(ai), 1 < i < m, independent 
and identically distributed exponential random variables then r^ a> _ >a y is the 

cumulative distribution of e~ £l ^ ■ e~ £2 ^ ■ ... ■ e~ £m ^ = e~ lm > a , where 7 m;Q 
is the gama random variable with the shape parameter m and the scale 
parameter a, i.e., it has the probability density function of the form a m / (m — 
1)! i"™" 1 e~ ax l( 0) oo)(a ; ); m particular, for a = 1 comp. Proposition 4 in Jurek 
(2004). 

In a proof of the above theorem the following identity, that might be also 
of an independent interest, is needed. 

Proposition 1. For different complex numbers Zj,j = 1,2, ...,n,n+ 1 we 
have equality: 

n -, n -. n -. 

Ej^— ( n rbhllT^— < 6 > 

i=l k=l\k^i t=l 

Equivalently, for any different complex numbers Zj,j = l,2,...,n 7 we have 
the identity 

n n 

e n ZJ ^ = !» f° rau zec > ( 7 ) 

• 1 ; -| j / - Zk Zi 
t=l k=l;kfi 

that can be regarded as a canonical decomposition of 1 as a sum of finite 
products of complex fractions. 
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Since the characteristic function of each v e ID is uniquely determined by 
the triple [a, R, M] from its Levy-Khintchine formula we will write formally 
that v = [a, R, M]; for details see the Section 2.1 below. 

If v = [a, R, M\ and jW(v) := [a^, i? {/3} , M^>] and 

bM /3 '■= / x ||x||~ 1_/3 M(rfx) G M d ( or a Banach space (8) 

then we have 
Q W = + l)" 1 ^ + & Mij8 ), = /?(2 + P)~ l R, 

Mlfi(A) = [ T t i /P M(A)dt = (3 [ Mis" 1 A) s^ds, for A e B , (9) 
Jo Jo 

where Bq stands for all Borel subsets of M. d \ {0} (or E \ {0}. 

With these notations Theorem 1 gives the description of random integrals 
(2) in terms of the corresponding triples. 

Theorem 2. For distinct positive reals (3i, 02, ■ f3 n , coefficients Cj tJl de- 
fined by (4), an infinitely divisible probability measure v = [a,R,M] and 
j{Pi,...,M^ = [ a {Pi,-,M ^ jj{Pi,...,M ^ m^'-'M] we have 

n „ n n t n p n 

« ,e w =«n^r+E^f- n ^=e^« m do) 

j=l P > + 1 j=l P > + 1 fc#,k=l Pk Pl 3=1 



n 



R { ^ = UJ^ R = J2 c ^R m , (ii) 

j=l 3 j=l 
r-1 pi n 

M&>~>M(A)= ... T Vft . 1/fSn M{A)dt l ...dt n = Y j C j>n M^\A) 
Jo Jo 1 "' n j=1 

(12) 

where bM,p j} a> ' 

are given in (8) and (9). 

Corollary 1. For different positive reals 0i, 02, Pn o,nd the constants C^ n 
given in (4) we have 

where for Cj >n < the corresponding convolution power means the reciprocal 
of the corresponding infinitely divisible Fourier transform. 

Here we have analogous formulae to that of Theorem 1, for the composi- 
tion of integral mappings with different collections of /3's. 
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Proposition 2. In the formula (2) in Theorem 1 we have: 
(a) For distinct positive reals p\, j3 n , a and integer m > 1, 

n 

m — times ^ ^ 

m—1 n 

- « _1 E ( E b d tn k) >{a, ay(t), < t < 1, (13) 

k=0 j=l 



fc + 1 — times 



where the coefficients d^ n are given by 



d«i = dg((3 1 ,...,(3 n] a^..^)= J] 



( A -ft a - & 

l-times ^rJ 

= C h „(/?!,.. .,/?„) (—^-5-)', Z = 0,1,2,.... (14) 
For different positive reals a and 7, and positive integers k,l > 1, 

k 

7 7-« L ^ (A;-s)! 

k — times I — times s — times 

k+l-l 

- - E W^^t^^W], /or < t < 1, (15) 

^ ^ r — times 



where to coefficients e r ^+iare given by 



s=l 



and /or 10GR and m — 1,2, (w) rn := w(w + l)...(u> + m — 1) denotes the 
Pochhammer symbol. 

In particular, we get an identity 

(—)*(—)' IE 7tH(— )*- E^.«(— r]" 1 < 16 > 

a — 7 7 — a L - iL ^ A; — s ! a 7 J 

s=l v ' r=l 

Remark 3. (i) Using the above formulae, and repeating elementary calcula- 
tions from proofs of Lemmas 2 and 3 one can get r& for the most general col- 
lections A = p 2 , ■■■fin, di,... a r , ..a r }, where p\, j3 n , a±, a r 

mi— times m r — times 
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are distinct positive reals. However, expressions are getting quite involved 
and are not presented here; cf. Section 2.7. 

(ii) For a finite set B of not necessary distinct complex numbers, let us 
introduce the function pi given by 



c^b, cGB 

Then for the set A from (i) we have 



n 



n< 



) mi for 1 < j < n 



«-)=n^- n 



for 1 < s < r. 



(18) 
(19) 



k=l ' l=is,l=l 

Furthermore, if we put K := n + mi + m 2 + ... + m r and define 



8j for 1 < j < n 

oti for n + 1 < j < n + mi 

a 2 for n + mi + 1 < j < n + mi + m 2 



a r for K — m r + 1 < j < K 
then A = {/3[, j3' 2 , j3' K } and hence we get that 



K 



ii a' — Ri 



for 1 < j < K, 



(20) 



which coincides with (18) and (19). The notation of pi ma y gi ve a more 
concise way of getting formulae for the general case of A; comp. above note 
(i). From the property of the constants C^ n in Theorem 1 we conclude that 



$^Pb(&) 



(21) 
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whenever B is a finite set of distinct numbers. 
2. Auxiliary results and proofs. 
2.1. Random integrals. 

Let us recall that for a probability Borel measures fi on IR d (or on E), its 
characteristic function ( Fourier transform) fi is defined as 

jj( y ) ■= f e i<y ' x> fi(dx), y E R d , (or y e E') 
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where < •, • > denotes the scalar product; (in case of Banach spaces, < 
•, • > is the bilinear form on E' x E). Further, the characteristic function 
of an infinitely divisible probability measure fi admits the following Levy- 
Khintchine representation 

jl(y) = e^ y \ y G M. d , and the Levy exponent $(y) = i < y,a > — 

\<y,Ry> + [ \e l<y ' x> -l-i<y,x> l B (x)]M(dx), (22) 

2 JR d \{0} 

where a is a shift vector, R is a covariance operator corresponding to the 
Gaussian part of jj, and M is a Levy spectral measure. Since there is a one- 
to-one correspondence between measures jj, G ID and triples a, R and M 
in its Levy-Khintchine formula (22) we will formally write fi = [a,R,M]; cf. 
Araujo-Gine (1980) or Parthasarathy (1967). 
Note that for s G K we have 

$(sy) = i < y, s(a + / x(l B (sx) - l B (x))M(dx)) > -^s 2 <y,Ry> 

JE\{0} z 

+ / [e i<y > z> -l-i<y,z>l B (z)]M(s- 1 dz) (23) 

JE\{0} 

Finally, let us recall that 

M is Levy spectral measure on R d iff / min(l, | |x| \ 2 )M{dx) < oo (24) 

For infinity divisibility on Banach spaces we refer to the monograph by 
Araujo-Gine (1980), Chapter 3, Section 6, p. 136. Let us emphasize here 
that the characterization (24) of Levy spectral measures is NOT true on in- 
finite dimensional Banach spaces ! However, it holds true on Hilbert spaces; 
cf. Parthasarathy (1967), Chapter VI. 

For this note it is important to have the following technical result: 

Lemma 1. If the random integral A = h(t)dY u (r(t)) exists then we have 
\og£(A)(y)= [ log C(Y^l))(h(s)y)d r (s)= [ <t>(h(s)y)dr(s), 

J(a,b] J(a,b] 

where y G M. d (or E' ) and $ is the Levy exponent of C{Y U {1)) = v. In 
particular, ifr is the cumulative probability distribution function of a random 
variable T concentrated of the interval (a, b] then log C(A)(y) = K[&(h(T)y)}. 
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The formula in Lemma 1 is a straightforward consequence of our definition 
(integration by parts) of the random integrals (**). The proof is analogous 
to that in Jurek-Vervaat(1983), Lemma 1.1 or Jurek-Mason (1993), Lemma 
3.6.4 or Jurek (1988), Lemma 2.2 (b). 

Remark 4. Note that for bounded intervals (a, b] C 1R + , continuous h and r 
of bounded variation, the integrals of the form A in Lemma 1 are well-defined. 



2.2. Proof of Proposition 1. 

Firstly, assuming the convention that Ylied Zi ~ ^ then the statement (6) 
is true for n = 1. 

Secondly, note that the last value z n+ i is not 'reached' by any other Zi, 
because 1 < i < n. It might be treated as a number independent of the index 
i. Hence, by an induction assumption, the equality (6) for different complex 
numbers z±, z 2 , .., z n _\ and the last values equal to z n and z n+ i, respectively, 
gives that 

n— 1 -. n—1 , n—1 -. 



. =i Zi z n \ k=ik _^.z k ZiJ k=i z k z n 



and 



n— 1 -. n—1 -. n—1 -. 

z^— ( n t ± t)=u^t- (26) 

Consequently, for different n + 1 complex numbers z±, Z2, z n , ^n+i we get 

n n 

^— ^( TT — ' 

Zi— 2 n + l y J. J. Z k —Zi 

1=1 k=l;ky^i 
n—l n—1 n 

=y— ^- f n — — -)+^- n — 

Zi— 2 n +l y J. J. 2fc— Zi 2 n — 2i y 2 n — 2 n + l J. J. 2 fc — 2 n 

i=l fc=l;fc^i fc=l;fc^n 

n— 1 n—1 n—1 



= n — U^^TT — 

Z / 2i— 2 n+ l 2 n — 2i ^ J. J. Z k — Zi J Z n —Z n + l J. J. 2 fc -2„ 

i=l k=l;k^i k=l 

n—1 n—1 n—1 

= Vf^ + ^)^— ( TT — V^- TT — 

Z_-/ \2i— 2 n + l Z n — Zi J Z n — Z n +\ y J. J. 2 fc — 2j y Z n — Z n +\ J. J. 2 fc — 2 n 



i=l k=l;k^i k=l 
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n— 1 n— 1 n— 1 n— 1 



2n— / / Zi—Z n + l ^ J. J. Z fc — Zj y 2 n — 2 n+ l / / Z n — Zj \_ J. J. Z k —Zi J 

i=l k=l;k^i i=l k=l;k=£i 

n-1 

+^hr n ^ ^ ( 2s ) and ( 2g ) 

fe=i 

n— 1 n— 1 n— 1 n— 1 

= ^^TT — —Y—( TT — H^^TT — 

Zn — Z„+1 A.A. Z n + \ z n — Zn+l Zi~ Z n ^ J. J. Zfe — Zj y Z n — Z n + l J. J. Z fc — Z n 



fc=l ^ n+1 i=l fc=l;fc^i k=l 



it ... 71—1 71—1 it 

TT — — TT — + ^^TT — = TT 



k' = -\ k' = -\ k = = 1 k^ = 1 

which proves the formula (6). 

By multiplying both sides of (6) by the product nr=i(- 2 * ~~ ^n+i) we arrive 
at identity (7) for arbitrary z = z n+ i different from all Zj, j — 1, 2, n. For 
z — zi, for some 1 < I < n, in (6) the term i = Z is equal 1 and all others are 
zero. This completes a proof of Proposition 1. 

Remark 5. (a) Note that it is enough to prove (7) for z — 0. This is so, 
because we may write (7) for a new constants z\ := z\ + z. 
(b) In the interpolation theory of functions for given set of points (xo,yo), 
(x 1 ,y 1 ), (x n , y n ) in the plane IR 2 , with distinct x , xi, x n , 

n n 

%k % 



p n ( x ):=J2yj II 

3=0 kjtj,k=0 



Xfc Xj 



is the unique Lagrange interpolating polynomial of degree less or equal n — 1 
and such that 

P n {xi)=yi for all % = 0, 1,2, 

cf. Kincaid and Cheney (1996), Chapter 6. Thus for the particular points 
(x , 1), (xi, 1), (x n , 1) in M 2 we get the line y = P„(x) = 1 as the Lagrange 
interpolating polynomial. However, our Proposition 1 is for complex coef- 
ficients with completely different proof and seems to be not known in the 
interpolation theory. 



2.3. Products of uniform distributions: part I. 

Here are some elementary identities concerning the products of powers of 
independent uniformly distributed random variables. The main objective is 
to express the cumulative distribution function (c.d.f.) or probability den- 
sity function (p.d.f.) of such products as a linear combinations (with not 
necessary positive coefficients) of other c.d.f. (or p.d.f.). 
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Lemma 2. Let Ui, 1 < i < n be i.i.d. uniformly distributed over the inter- 
val (0,1], a,i > and let /{ ai)tt2 ,..., a „} and F {ai)a2) ... )Qn} denote tee probability 

density and cumulative distribution function of U^ ai ■ U 2 ^ a2 ■ ... ■ U}J an , re- 
spectively. Then 

( a ) f{ai,a2,...,a n }(.%n) = 

J Xn J X n — i J X n —2 

... I a 2 xf- a ^ 1 j a 1 x c f 1 - a2 ~ 1 dx 1 dx 2 ...dx n _ 2 dx n _ 1 , < rr n < 1. (27) 
//ai = a 2 = ... = a m = a then 

fi^ayix) = ax*- 1 t^^l for 0<x<l, (28) 

m — times 

m—1 / i xj 

F {a ,...,a } (s) = s a J2 jr M < s < 1, (29) 

m— times ^ ^ 

and F{q; . a}(s) = 1 for s > 1 and zero /or s < 0. 

m-times 

^ If all positive reals cti (i = 1,2, ...,n) are distinct and 

n n \ 

C j>n (ai, ...,a n ) = C jtn : = ] [ -j— and Cj, a := ]] — -(30) 

then 

n n 
f{ ai ,a 2 ,...,a n }(x) = C 0,n CtjX*^ 1 = «!...«„ Cj,^" 1 , < X < 1, (31) 

i=i j=i 

and 

i ? {ai ) a 2l ... 1 a„}(s) = 5^ Cj,n s a '' , /or < s < 1, wnere C j>n = 1. (32) 
i=i j=i 

Proof. For positive and independent rv X and Z with p.d.f. fx and /z, 
respectively we have that X ■ Z has the p.d.f. 

/*.z(*)= / fx(-)-fz(x)dx. (33) 
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Since f{ a y(x) = ax a 1 l( 0jl )(x) is the p.d.f. of U 1 ^ therefore from (33) we 
get 

f{a,P}( z ) = fw/°-uw( z ) = P / f{a}{-)^~ 2 dx = f{f j} (z) / f {a} (t) t~ p dt. 

J z x J z 

(34) 

Hence for a± and a 2 we conclude that 



f{ai,aa}M = a 2 X% > 1 / ftl^ 1 ° 2 1 dx x 

J X 2 



(35) 



which is indeed of the form (27) for n = 2. 

Assume, by the induction argument, that the formula (27) holds true for 
n. Then using (33) and (34) we obtain 

f{ai,a2,-..,a n ,a n+ i}{ x n+l) ~ fa n +i{ x n+l) I f{a 1 ,a2,...,a n }(x n )x n + dx n — 

Jx n+1 

Qn+i-1 f ( a„-l / a„-i-a„-l /" a„-2-a„_i-l 

«n+l x n +l / \ a n-L n / "ra-l x n-l / "re-2 x n-2 

Jx„ +1 V A„ J 'x„_i 

J X3 J X2 

^n+l^n+l / X n + / ^n— 1 X n— 1 / ^n— 2 -^n— 2 

Jx n+ i J x n J X n -l 

... / a^a^ 2- " 3 " 1 / «i re" 1- " 2-1 cfai drr 2 ••• dx n -i dx n , 

J XZ J X2 

which is the equality (27) for n+ 1. Thus the proof of the part (a) is complete. 

Taking in (27), a\ — a 2 — ■■■ — oc n — a and performing the successive 
integrations, we get the formula (28). (Or simply prove (28) by the induction 
argument utilizing (33)). Integrating p.d.f. (28) we get c.d.f. (29). (Recall 
also here Remark 6 and the formula (39) from below.) Thus the part (b) is 
proved. 

In the part(c), formulae (31) and (32) are obvious for n=l. Assume that 
(31) holds true for n. First, from Proposition 1 formula (6) we infer that for 
1 < j < n we get 

n 

Cj,n{0tn+1 — Oij) 1 = Cj >n+ i, (cKj — a n +l) lc j,n = Cn+l,n+l ■ (36) 

J'=l 
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Then from (34), (31) and (36) and again (6) from Proposition 1 we get 

f{a 1 ,a2,...,a„,a„ + i}(z) =a n+l x ™ +1 / /{ai, ai, ...,a„} {t) t n+1 dt 

J x 



n „i 

= ai...a n+ i °j,n / t^-^-'dt 

n 1 



V c J>+ i^ 1 + ai...a n+ i ( V c j)n )x a ^ 



n n+1 

= a 1 ...a n+1 (^2 c jtn+l x a i- l +c n+1)n+l x an + 1 ' 1 ) = ai...a n+ i c jjn+1 x aj ~\ 

3=1 3=1 

which completes a proof of (31). The expression (32) for c.d.f. is immediate 
consequence of (31). This completes a proof of Lemma 2. 

2.4. Proof of Theorem 1. 



In view of Lemma 1, to prove formula (2) it is necessary and sufficient to 
show the equality 

\og({J^M{v))){y) = E[\ogu(Ul /h ■ Ul /h • ... • U^y)} 

= f log u(ty)dr {Pl _ M (t), yeEf. (37) 
J o 

Of course, (37) holds for n — 1. Assume it is true for n—1. Then from Lemma 
1 and the definition of the mapping (given before Theorem 1) we 

get 

log(J { ^'-'^- lM (u)j(y) = f \og(J { ^-^\v))(sy)ds^ = 

Jo 

f\[\ogv{ul lh .ul /h -..,U^r 1 sy)\ds^ = E[\ogv(ul /h -U 1 2 /h -..,U^y)} 
Jo 

= / logv{sy)dr { p u ...,{ in _ M (t), (38) 
Jo 

which completes proof of (37) and consequently the formula (4). Explicit 
expressions for time changes 7"{/3i,...,/3 n _i,/3 n }(i) are given in Lemma 2, part (c). 
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2.5. Proof of Theorem 2. 



First, putting $(y) = log z>(y) ( the Levy exponents of v) into (37) we get 

i < y,a {l3l '-' M > ~ < y,R {(3l '-' M y > 

+ / [ e i<v,*> _ i _ i < y> x > l B (x)]M {/3l -- /3 " } (rfa;) 

</R d \{0} 

= E[$(^ 1//3l -f/ 2 1/& -...-^/^)] 



Since for < s < 1, we have that 1b(sx) — 1b(x) = l{i<||rr||<s-i}(^) therefore 
from the above and (23) we get 



a {/3l,-,/3n} 



= E[U 1 1 ^-U 1 2 ^:...U 1 n ^(a+[ xM(dx))} 



IlTX^ a+ / / srrM(^)dr {/3li ... i/M (s) (by (4)) 

J=l 1 + Pj ./0 ^K||x||<s-i 

n^T^^-A)]^ / / sxM(dx)s^- 1 ds (by (8)) 

,=i 1 + Pi 7= i -/o 7i<||x||< s -i 



=ni^-+E^i( n ^ (7) z = -i) 



j = l J fc^j,fe=l 

j= i + Pj j= i k &,k=i Pk p i j=i 

which proves the formula for the shift vector. 

Similarly, for the Gaussian part, using again the identity (7) (with Zj = (3j 
and z = —2) we get 

R Wl,..,M = E[([/ l/A . fjl/fh . . Jjl/Pn^ R = JJ R 

3=1 ' j + 2 

n n n n n . rt n n q n 

=n^(E n n i^)** 

i = l ^ + 2 J=l fc ? 4i ) fc=l « |=1 fc # i,fc=l P fc « J=l 

which gives the formula for Gaussian covariance. 
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Finally for the Levy spectral measure using (4) and (9) we have 
M&>->M(A) = E[M((U} /pl • U] IP2 ■ ... ■ U l J^)- l A)} 

= / T s M(A)dr {Pl ,..., M (s)^/3 1 /3 2 .../3 n J2 c J,n / M(s' 1 A)s^- 1 ds 
Jo j=1 Jo 

n n p. „i n 

= BE ir\n M(s-^A)ds = J2c„ n M^(A), 



which completes the proof of Theorem 2. 

2.6. Products of uniform distributions; part II. 

For the proofs below note that for a ^ we have the identity 

/VM-log^-W=<^[l-/£fc^], .>0, (39) 
^ S a k=0 

because both sides vanish at s = 1 and have the same first derivative. 

Remark 6. Note that for a > 0, the expression on the right hand side of 
(39) in the square bracket, coincides on (0, 1] with the c.d.f. -F{ a ,...,a} from 
(29); recall also Remark 2. Below we will use the same notation F{ 0j ... ja } for 
a < but will not treat them as the c.d.f. of products of negative powers 
of independent uniform r.v's. Nevertheless first derivative of F{ a ^ is still 
f{a,...,a} given by (28). 

Lemma 3. (i) For distinct positive reals (3i, (3 2 , (3 n , a let us define 

l~times 

n ^(^t)'-^a m^-;. '-m,2,.... (40) 

k^j ,k=i 1 K • 3 13 • 3 

Then for n > 1 and m — 1, 2... we get 

n m—1 n 

f{^,...,0 n ,a,...,a}(x) = J2 d fn Yl (J2^ d l>~ k) )f{^,-,a}(x), 

3= i &=o ?'=i l ! " Y" 

** J fe + l— times 

(41) 
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and hence 



m—l n 



y 1=1 fc=o 7=1 . ', v . . 

m — times J J k + l — times 

n n 

= E d fn F{^}{x)+F { a 1 a l ^_ 1 a}{s)- [E 4> F { a - a - fy } F w( s ) ■ 

3 ^ m — times 3 ^- 

m — times 

(ii) For k, I > 1 and positive j ^ a we have 
f{a^^,^-^J}( z ) = (- - (- ^ 

k — times I — times 
k /,\ k + l — 1 



s — times r — times 



where coefficients e r ^+i ar & given by 

rAk / \ /,\ rAfc 



^-L(-D (P _ ()l(t _ ()l -L(- 1 ) U - 1 A '-i J' 

/or l<r</c + / — 1 and for wGl and m — 1,2, (w) m :— w(w + l)...(w + 
m — l) denotes the Pochhammer symbol. 

Proof. Using (33), (31), (39), (40) and then (28) we get 

/{ft,... A.ai^^}^) = y f{/3 1 ,...,/3 n }(j)t~ 1 f{a,...,a}(t)dt 

m — times 

n /m n pi 

= (^3T)T E r - 'V'-' '/. ^-'(-logO™- 1 * 

- V" >) f ^-i v- Cj.n gj {{fa - a)\ogx) k ^ 

- 2^ d j,n tm X. {a _ « r jfci ) 

j=l j=l K ' ■>> k=0 



n m—l , , ■ i. a 

j=l k=0 ' j'=l V 



" m—l / , \k 11 

i=i fc=o * v : 'j=i 

n m—l n 

= E 4? Aftl - « _1 E (E & d tn k) ) ha, a, a } {x) 



(m—k) 



i— 1 fc_0 3— 1 fe+ i_ii mes 
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where the probability density functions f^.y and f{ a , a ,...,a} are given by (31) 
and (28) respectively. Thus we proved the formula (41). 

From the above immediately follows the c.d.f. i^^,^,...,^^,...^} i n the 
first form in Lemma 3 (i). For the second expression note that 



m— 1 n 



fc=0 3=1 




j=l k=o ^ 1=0 

n m—1 , , \t m—1 

3=1 1=0 ' k=l * i 

n m-1 , , \i m-1 , , \ ; 

V-^ V- (-a logs)' „ a Nm _; v „ _ Q ( z abgsy 

j=l Z=0 ^ «=0 

-2^^>(^r^) n F{a} 

j=l 10 1=0 

n 

= i \a,a,...,q }(s) - ^ C ^WZrj ^ F {&}( s ) F {a - ft, ...,a - Z^ 5 )' 
^T^ZT i =1 J v v ' 

m — times 

which completes the proof of the part (i). 

Before the proof of part (ii), for an ease of notation, let us recall here the 
Pochhammer symbol 

for a e Rand n = 1,2, ... (a)„ := a(a+l)(a+2)...(a+n-l), (a) := 1. (42) 

For fc > 1, Z > 1 and 7^a and 7, ct > using (33) and (39)we have 

k — times I — times k — times I — times 

^—t " a - 1 f x^- 1 (- log z + log x) fc_1 (- log x)'- 1 ^ 



(fc - 1)!(Z - 1)! 

E(-!) s (* 7 1 ) { ~ logz)k ~ 1 ~ s J! x ' r ~ " 1 (- ]o & x ) 



(k-m-iy. o 
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(* " !)!(' - 1)! H V 3 J (1- Q ')* + ' 

•['-^ H, -; 1W ] 



i=0 



fe-1 



\fc-l-s 



= °* (— )' E ^~T7 — ~ — r" „ ^Ml-^ - a 7 - «,(*)] 

7 _a ^ s! («-7) s (A;-l-s)! V 7 a ) y n 

l-\-s — times 

k — 1 

= (77V E % % 3>(^)[ 1 - ^ - a, 7 - 

k-s-times l + s-times 

= (^)'(^)' E .... y WH-iy - a. ....T - qMl 

s times ;_[_fc_ s _^ mes 

(43) 

For further simplification of (43) note that the function (28) and (29) are 
well defined for a, b G H. \ {0}. ( See also Remark 6.) Hence, if a + b ^ and 
p, q > 1 then 

p-times q-tim.es 3 p+j-times 



j— times 



U - p) 

Using the above with p = s, q = I + k — s, a = a and b = 7 — a we get 
f{a, ...,a}(z)[l - -F{7 - a , ...,7 - a}(^)] 



Z + fe — s — times 



^ l+k-l >y — a 

s — times S r — times 

Hence and from (43) we get 

k 

h^^^Mz) = (— ) (— ) ^ -) 



k — times I — times 



[%, .... ,9 w - (-^r 'e' #%t( V r %-iz>m i- < 45 > 
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Changing the order of summation in the second summand in (45) we get 

s=l K y r=s x 7 Tt 

r — times 

r=\ S =i v 7 ' v >' ' y 

r — times 

+ E 1 [E(-i)'fe f fe f ](^)'/^w. 



r=fc+l s 1 

Defining the coefficients 

= f EI =1 (-l)^f^ for 1 < r < k 
^ 1 ELi("l)fe {fet farfc + l<r<fc + Z-l, 

finally give us that 



k — times I — times 



k+l-l 



IE ra(V r/ <t-22> (2) - E ^(^r^wi. 



which completes a proof of (ii) in Lemma 3. 
2.7. A note on a more general case. 



For distinct f3i, /3 n , a, 7 using again (33) and Lemma 3 we obtain 

f 1 z 1 

/{ / 9 1 ,/3 2 ,..,^ n ,a^^,^V^j}(^) = y /{ / 9i,..., / 9„,o^^}(-)-/{qv^}(^)^ 



m — times I — times m — times 

m—1 n 



3 ^~ I — times k j ' 1 fc+l-hmes I — times 

(46) 

where p.d.f. /{^^^(z) and /{qi^^qv-^J}^) are given in Lemma 

I — times k-\-l — times I — times 

3. Thus, in principle, for an arbitrary finite set A of positive reals we have 
a formula for /a(z). However, our expressions become quite complicated; 
comp. Remark 3. 
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3. Concluding remarks. 

(a) The formulae for the time scale change r in Theorem 1, Proposition 
2 and in the section 2.7 are quite involved and not very transparent. Maybe 
using the notation pb given in Remark 3, formula (17), would clarify the 
expositions. 

(b) The ideas of the proof of Theorem 1 seems to be applicable to com- 
positions of mappings I^b]i wnere the time scale change r is a c.d.f. Namely, 
from Lemma 1, if r\ and r 2 are c.d.f. of rv's X±, X 2 concentrated on intervals 
(ai,6i] and (02,^2], respectively, then 

log (^(CSiWKy) = Epogi>(M*i)M*2)y)] 

hg Hsy)dr 3 (s) = \og{ll2 b3] ^)), (47) 

(03,63] 

where (a 3 , 63] = r 1 ((ai,bi\) ■ r 2 ((a 2 ,b 2 ]) an d r 3 is the cumulative distribution 
function of rv hi(Xi)h 2 (X 2 ). However, if hi(xi) ■ h 2 (x 2 ) = hs(g(xi,x 2 )) for 
some measurable functions g : R 2 — > R and /i 3 : R + — >■ R then 

rhi,ri T h 2 ,r 2 _ T h 3 ,fi / \ 
i (oi,6i] ° I (a 2 ,b 2 ] ~ J (a 3 ,63]V^' 

where f 4 is the c.d.f. of rv g(Xi,X 2 ) and o denotes the composition of 
the random integral mappings (provided the mappings in questions are well 
defined, i.e., appropriate domains and ranges coincide). 
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